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Abstract. - We derive the low temperature properties of the molecular magnet {Mo72Fe3o}, 
where 30 Fe^+ paramagnetic ions occupy the sites of an icosidodecahedron and interact via 
isotropic nearest-neighbour antiferromagnetic Heisenberg exchange. The key idea of our model 
(J.S. & M.L.) is that the low-lying excitations form a sequence of "rotational bands", i.e., for 
each such band the excitation energies depend quadratically on the total spin quantum number. 
For temperatures below 50 mK we predict that the magnetisation is described by a staircase 
with 75 equidistant steps as the magnetic field is increased up to a critical value and saturated 
for higher fields. For higher temperatures thermal broadening effects wash out the staircase and 
yield a linear ramp below the critical field, and this has been confirmed by our measurements 
(R.M.). We demonstrate that the lowest two rotational bands are separated by an energy gap 
of 0.7 meV, and this could be tested by EPR and inelastic neutron scattering measurements. 
We also predict the occurrence of resonances at temperatures below 0.1 K in the proton NMR 
spin-lattice relaxation rate associated with level crossings. As rotational bands characterize 
the spectra of many magnetic molecules our method opens a new road towards a description 
of their low-temperature behaviour which is not otherwise accessible. 



Introduction. - A new class of magnetic compounds known as molecular magnets [y is 
attracting much attention. These compounds can be synthesized as single crystals of iden- 
tical molecular units, each containing several paramagnetic ions that mutually interact via 
Heisenberg exchange. The intermolecular magnetic interactions are in the great majority of 
cases utterly negligible as compared to intramolecular magnetic interactions. Measurements 
of the magnetic properties therefore reflect those of the common, individual molecular unit. 
Molecular magnets such as {Mni2} and {Fes} have been the focal point for intense study 
of subjects of broad scientific importance, such as quantum tunneling of magnetisation and 
quantum coherence j2| . However, it is doubtful |^ that larger symmetric arrays of paramag- 
netic ions can be accommodated in such polynuclear coordination complexes utilizing simple 
organic bridging ligands. Very recently the first examples of a new paradigm of molecular 
magnets, based on so-called Keplerate structures, have been synthesized Q, and these offer 
numerous avenues for obtaining truly giant, highly symmetric arrays of paramagnetic ions. 
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The archetype of this new class is referred to as {Mo72Fe3o} [Q. Embedded within a (diamag- 
netic) host molecule (mol. wt. 18,649), 30 Fe^+ paramagnetic ions (spins s = 5/2) occupy 
the sites of an icosidodecahedron and interact via isotropic, nearest-neighbour antiferromag- 
netic exchange. This Keplerate and its interlinked derivatives pose a major theoretical 
challenge. The dimension of the Hilbert space for {Mo72Fe3o} is a staggering 6"^°, precluding 
the calculation of the energy eigenvalues on any computer. 

In this Letter we show that despite this immense obstacle the major low-temperature 
properties of {Mo72Fe3o} can be calculated. Rather than pursue the futile task of diagonalizing 
the Heisenberg Hamiltonian [see (Q)], we adopt an approximate, diagonalizable Hamilton 
operator [see d)] which properly incorporates the generic result that the low-lying magnetic 
energy levels of a wide class of molecular magnets are arranged as parallel rotational bands. 
These rotational bands reflect the underlying sublattice structure of the spin array |^,|^ that 
has been found for the corresponding classical system [^ . But in contrast to the classical model 
and approximations like high-temperature series expansions the present approach enables us 
to access the properties of giant magnetic molecules for very low temperatures. We are also 
able to show that the predictions of the classical model extend to temperatures as low as 
100 mK with very good accuracy for the special case of {Mo72Fe3o}. In particular, we present 
our theoretical predictions for the temperature and field dependence of the magnetisation, 
and for the resonances of the spin-lattice relaxation rate of {Mo72Fe3o}. We also show that 
our experimental findings directly reflect the existence of rotational bands. 

Rotational band Hamiltonian. - The Hamiltonian for the isotropic Heisenberg model 
including the interaction with an external magnetic field B reads 

H ^ -2J Mu)-l{v)+gfiBBS, , (1) 

{u,v) 

where J is the exchange constant with units of energy, and J < corresponds to antiferro- 
magnetic coupling, g is the spectroscopic splitting factor, and /is is the Bohr magneton. The 
vector operators s (u) are the spin operators (in units of h) of the individual paramagnetic ions 
with spin s. The sum in (^ runs over all distinct nearest-neighbour pairs (u, v) of spins of a 
single molecule at positions u and v. Since the Hamiltonian commutes with the operators 
and Sz of the total spin, the eigenstates of H may be classified using the quantum numbers 
S and M. 

In order to develop an approximate quantum model of {Mo72Fe3o}, we first exploit the 
fact that the set of minimal energies for each S forms a rotational band, 

Es,mu. « -J [DiN, s)/N] S{S +1) + Ea. (2) 

This has been noted on several occasions for ring structures with an even number N of 
sites Elsewhere we reported that all finite Heisenberg systems with isotropic and ho- 

mogeneous antiferromagnetic exchange (including rings for both even and odd N , tetrahedron, 
cube, octahedron, icosahedron, triangular prism, and axially truncated icosahedron) exhibit 
a rotational band [^ . For high-symmetry systems we provided an expression for the classical 
limit D{N, oo), henceforth denoted by D, which yields D = A ior rings with even N, see 
also [|l^, and 13 = 6 for the icosidodecahedron, cube, and octahedron. Our investigations [^ 
have shown that the numerical value of D(N, s) for any finite s is always a little larger than 
D. This difference can be taken as a measure of quantum spin effects; consistently it is largest 
for 3 = 1/2. 
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The constant offset Ea in (|^) is to be selected so that the highest level of the rotational 
band, which occurs for 5* — Ns, agrees with the largest energy eigenvalue of (0). The corre- 
sponding eigenvector is the ground state eigenvector for the counterpart ferromagnetic system. 
Thus the largest energy eigenvalue serves as an anchor in all approximations of the rotational 
band Summarizing, for the molecular magnet {Mo72Fe3o} the lowest rotational band is 
given by 
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y-sis- 
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Fig. 1 - Low-lying energy eigenvalues using the approximate Hamiltonian with B = 0. The levels 
of the lowest band are given by (b|) . 



The rotational band (||) by itself is insufficient for calculating observables at low temper- 
atures; it is essential to include the other low-lying excitation energies. This can be achieved 
via the introduction of an effective Hamiltonian whose form is set by incorporating the known 
symmetry of the spin array (icosidodecahedron) , which is similar to that of the triangular 
lattice. Therefore, the low-lying spectrum can be understood as originating from three in- 
teracting sublattice spins [Q. These findings are supported by the fact that also the exact 
classical ground state of the icosidodecahedron for J < is describable in terms of three 
sublattice spin vectors Sa, Sb, and Sc with Sa ~ Sb ^ Sc = 25, and with relative angles of 
120° [p|. We thus adopt as our effective Hamiltonian, replacing the field-free term of (|l|), see 
also [HT7|,|ll|l, 



^1 



cff 



J 



S'4 



Si 



si 



(4) 



The three sublattice spin quantum numbers, Sa, Sb, Sc, can assume values 0, 1, ... , 25. The 
sublattice spin operators mutually commute and they also commute with Hl^. Thus the 
eigenvalues of H^^ are given in terms of the quantum numbers S, Sa, Sb, and Sc- The 
lowest rotational band, (||), arises upon choosing Sa — Sb — Sc ~ 25 and allowing S to 
extend from to 75. The next higher rotational band, see Fig. |l|, is obtained for the choice 
Sa = Sb = 25, Sc = 24 and its permutations. Note that these two bands are separated by an 
energy gap A — lOJ. Continuing this process leads to a sequence of parabolic bands. This is 
rather realistic for the second band, and indeed observed in many finite Heisenberg systems, 
compare Figs. 1, 3, and 4 in Q. 
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Experimental implications. - The theoretical result for A4/{gfiB) versus B at T — K, 
as obtained using (Q) and a Zeeman term, is a staircase with 75 steps of unit height which 
terminates at the critical field Be = 30\J\/ (g^B) = 17.7 T. The values of J and g have 
been determined by high-temperature, low-field magnetic susceptibility measurements p^ : 
J/ks — —0.783 K, and g — 1.974. For B > Be all spins are aligned parallel and the total 
moment is given by = Tbg^B- 




1 2 3 5 10 15 20 25 30 

B (T) B (T) 

Fig. 2 - L.h.s.: The magnetisation according to is shown for T = 1 niK (staircase) and 100 mK 
(thick solid curve) as well as for the corresponding classical Heisenberg model at T = (dashed curve). 
R.h.s.: The thin curve displays the magnetisation following from (^) for T = 4 K. Experimental data 
using a pulsed field are given by the solid squares, their size reflects an uncertainty of ±0.5 T for the 
data. The thick curve is gives the result of the improved approximation (|^, also at 4 K, and closely 
reproduces the measured values. 



In the left panel of Fig. || we provide a graph of M/{g^B) versus B as predicted by the 
quantum model for T = 1 mK (staircase) and 100 mK (thick solid curve), as well as the 
rigorous result for the corresponding classical Heisenberg model at T = 0, a strictly linear 
dependence of M on i? up to Be (dashed curve). In fact, for T > 50 mK thermal broadening 
effects smear out the staircase of the quantum model and, above approximately B = 0.5T, A4 
increases linearly with B, almost parallel to but somewhat below the classical T = result. It 
would be extremely interesting to put this prediction of staircase-like behavior below 50 mK 
to a careful experimental test. 

What is currently available are experimental data for A4 versus B on {Mo72Fe3o}, see the 
right panel of Fig. ^, as obtained using a pulsed magnetic field. Due to the rapid introduction 
of the magnetic field, increasing from to 60 Tesla in approximately 6 ms, we estimate that 
the effective temperature of the spins was 4 K, even though the nominal cryostat temperature 
was 0.46 K. Nevertheless, the very good agreement between theory and experiment confirms 
the underlying picture of rotational bands on which the model is based. For a non-parabolic 
dependence of the energy eigenvalues on S the magnetisation curve would show unequal 
steps at T = as well as nonlinear behaviour at higher temperatures. As remarked above, 
the classical coefficient D = D{N, oo) always underestimates the true coefficient D{N, s) by 
a few percent Of course, there is no simple method for establishing the correct value of 
1^(30, 5/2), so it is of interest to estimate its value by using the currently available experimental 
data for A4 versus B. We attempt to improve approximation (^ quantitatively while ensuring 
that the rotational band structure is not altered. We adjust D{N, s) so that the resulting 
magnetisation curve (thick curve in Fig. 0, r.h.s.) provides an optimal fit to the measured 
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data. This is achieved by taking D{N,s) 
alternate effective Hamiltonian is given by 



6.23, which is very close to D = 6. Thus an 



Hf = -J 



N 



s'-i s\ + sl + s 



c 



(5) 



with 7 = 1.07 in order to maintain the correct value of the largest energy eigenvalue. Using 
(||) we find that Be ~ 18.4 T. Finally, using this tentative estimate for 13(30, 5/2) enables us 
to predict the ground state energy of {Mo72Fe3o} as Eq/Ub ~ —339 K. Taking into account 
the error bars of the measured data we estimate that the resulting uncertainty of 13(30, 5/2), 
7, and Eq is less than 3 %. We are currently attempting to calculate 13(30,5/2) from first 
principles using DMRG techniques [nsl. 
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Fig. 3 - L.h.s.: Magnetisation vs. temperature using the four lowest bands of the improved ap- 
proximate Hamiltonian (^ with B — 7.0 T (solid line); experimental data (dashed line). The near 
constancy of Ai with T provides strong evidence that the lowest rotational bands are indeed parabolic. 
R.h.s.: Low- lying energy eigenvalues in the range shown using Hf with B = 7 T. 



We now compare our theoretical results using (|^) and a Zeeman term to our magnetisation 
measurements at fixed magnetic fields. For both measured field strengths B — 5.5 T (not 
shown) and B = 7.0 T (Fig. |, l.h.s.) we find good agreement to the experimental data. 
The seemingly unexciting result, that both the experimental and theoretical magnetisation M 
are virtually constant over a wide temperature range, actually provides important, indirect 
confirmation of the existence of a sequence of rotational bands. This can well be understood 
with the help of the spectrum of Hf shown in Fig. ^ (r.h.s.) for a finite magnetic field. The 
levels of the lowest parabola in Fig. || (r.h.s.) originate from the M — —S levels of the lowest 
rotational band for i? = (Fig. |l|). Those of the second parabola in Fig. || (r.h.s.) originate 
from the M — —S + 1 levels of the lowest rotational band as well as the AI — ~S levels of 
the second rotational band. Note that the levels of the parabolas are distributed in a very 
symmetrical manner about their common minimizing value of S for the given B. Therefore 
these levels are populated symmetrically and this leads to a magnetisation which is nearly 
temperature independent. A slight temperature dependence for T>2 K (see inset of Fig. ^ 
l.h.s.) stems from small differences in the degeneracies of the individual levels as well as from 
the fact that the approximate spectrum deviates from the exact one for higher excitation 
energies. Nevertheless, the near constancy of M with T is therefore also strong evidence that 
not only the lowest rotational band (||), but also the second is parabolic, for otherwise the 
symmetry would be broken and M. would vary significantly with T. 
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Fig. 4 - L.h.s.:Energy levels of the lowest rotational band as functions of magnetic field. R.h.s. Low- 
temperature transitions, labelled 1 and 2, which will give rise to peaks in an EPR spectrum involving 
the lowest Zeeman levels of the two lowest rotational bands, separated by an energy gap Eg- The 
gap size Eg is independent of S, because the rotational bands are parallel. 



The existence of rotational bands implies also the occurrence of resonances of the spin- 
lattice relaxation rate at low temperatures iQ. Proton NMR probes whether two levels 
are separated by fko^ where is the proton Larmor frequency. Since this energy is small 
compared to the spacings between successive levels of the rotational band in the absence of a 
magnetic field, resonances of the relaxation rate could be seen close to level crossings induced 
by the static magnetic field In the absence of a crystal field the level crossings will be 
observable at the magnetic fields highlighted in the left panel Fig. ^. The zero-field level 
spacings of the lowest rotational band are As = 2(5' -I- 1) • 0.16 K, and thus it will be necessary 
to go to very low temperatures in order to observe distinct resonance peaks. 

We now discuss some of the other consequences of the existence of a sequence of rotational 
bands. Using Eqs. (^) or (^) we estimate that the gap between the lowest two rotational bands 
is A « 8 K. Inelastic neutron scattering and EPR techniques should provide useful tests of 
this prediction. In the standard EPR setup, working with a fixed resonance frequency, vq, and 
an applied field of variable strength B, because of the selection rule AS = one can focus 
on the Zeeman levels of the lowest two rotational bands for a fixed value of S. In the right 
panel of Fig. || we illustrate two distinct cases where the resonance condition is met for a given 
value of B. The arrow marked 1 applies for a case where the inequality vq < < 2vg is 
satisfied, where vq ~ Eq/H is the gap frequency associated with successive rotational bands. 
When the resonance condition hv^ — Eq + gfJ-sB is met (independent of S) the additional 
selection rule AM = ±1 allows for a transition from the M — —S level of the lowest rotational 
band to the M = — S* + 1 level of the first excited rotational band. For the arrow marked 
2, which illustrates the case vq < vq, the resonance condition reads hvQ = Eq — g^sB, and 
it is associated with a transition from the M = —S + 1 level of the lowest rotational band 
to the M — —S level of the first excited rotational band. The intensity of the peak will be 
temperature dependent, proportional to the Boltzmann factor exp{—f3g^BB). 

Anticipating that the Debye temperature of {Mo72Fe3o} is of order 200-300 K, we believe 
that specific heat measurements in the temperature range below 0.75 K can serve as a useful 
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probe of the levels of the lowest rotational band. A comparison between our theory predictions 
and specific heat data will be given elsewhere jl5j . 

Summary. - In summary, in this Letter we have presented an approximate quantum 
model enabling us to determine the low-temperature properties of the molecular magnet 
{Mo72Fe3o}. Several of our results have been confirmed by our measurements and others 
remain to be tested. The form of the effective Hamiltonian, Eqs. or (js]), is set by incorpo- 
rating the symmetry properties of the spin array which leads to an identification of interacting 
sublattice spin vectors. Some features of the model (e.g., parallel rotational bands also for 
high- lying states) are oversimplified, however, these do not play any role on results shown here. 
The present method, which is based on the notion of parallel rotational excitation bands, of- 
fers an insightful and quantitatively useful platform as an alternative to the insurmountable 
difficulties in treating the exact quantum Heisenberg model. Adaptation of this methodology 
to still larger, high-symmetry magnetic molecules that will be synthesized in the future can 
be expected to provide similar benefits. 
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